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TRANSMISSION LINE DESIGN FOR A POWER DISTRIBUTION SYSTEM 
AT 20 kHz FOR AIRCRAFT 

Leon W. Zelby, John B. Mathes, and John W. Shawver 
University of Oklahoma 
Norman, Oklahoma 73019 


1.0 INTRODUCTION 

The purpose of this work is to design a transmission line for a power 
distribution system in aircraft. The line is to operate at 20 kHz, have very 
low inductance and, consequently, a low characteristic impedance; it must with- 
stand 440 V at an altitude of 15 km (50 000 ft) and 1 kV at 300 km (about 
100 000 ft), and must be capable of carrying 100 A. The connectors and sup- 
ports must have a minimum creepage path of 0.635 cm (1/4 in.). 

Additional considerations of the line design include connectors for par- 
allel loads, number of such connectors, mechanical strength and heat dissipa- 
tion, dc resistance, and weight. 

The fundamental problem in determining the inductance per unit length, and 
the characteristic impedance of the line is the calculation of inductance or 
capacitance because LCv^ = l, with L and C representing the inductance and 
capacitance per unit length, respectively; and v representing the velocity of 
wave propagation along the line (in the case of the fundamental mode, the TEM 
mode, the velocity is governed by the properties of the material between the 
wires). At a frequency of 20 kHz, the skin depth for copper (resistivity 
1.7xl0 -8 fl-m conductivity = 5.88xl0 7 mho/m) is 0.46 mm n 0.5 mm, and the wave- 
length in free space, is 15.0 km and about 9.0 km in many commercially used 
insulating materials (relative dielectric constant on the order of 3). Thus, 
the total length of the line is on the order of 0.02 to 0.033 wavelength which 
justifies the quasi-static approximations (refs. 1 to 5). The close spacing 
among the individual conductors strengthen the validity of the relation 
LCV 2 = 1 and of the "quasi-TEM" mode dominance. [Strictly speaking, the TEM 
mode in this case is not absolutely pure because of two factors: there are two 

dielectrics between the wires, the insulation and air; and the wires are lossy. 
As a consequence, there will be an axial component of the electric field, but 
it is so small that the assumption about TEM is more than satisfactory for 
practical purposes. Furthermore, with reasonable insulation the shunt conduct- 
ance (between the wires) will be on order of 10~^ mho/m (ref. 6). As a con- 
sequence, the term RG/u^ lc will be on the order of 10 -14 « 1, so that 
the propagation constant will very nearly be equal to «VLC, so that the 
relation LCv 2 = 1 will be satisfied in any practical system.] 

The SI system of units will be used throughout this report. 


2.0 DETERMINATION OF PARAMETERS 

2.1 Introduction 

In order to design a suitable transmission line, several parameters need 
to be determined for different configurations. The parameters are L, the 
inductance per unit length; C, the capacitance per unit length; Z 0 , the 



characteristic impedance; and R, the resistance per unit length. (G, the 
shunt conductance per unit length is too small to be included.) L and C are 
related as indicated above through the relation LCv 2 = 1; also, for a lossless 
line, Z 0 = (L/C) 1 / 2 ft. At a frequency of 20 kHz, the skin depth in copper 
is 0.467 mm, with the current distributed in a circular cylindrical conductor 
as shown in figure 1. The distribution is governed by the relation (ref. 7) 



where Ber(x) + jBei(x) = J 0 (3 ~^ 2x )» with J o( x ) Bessel function of the 

first kind and zeroth order, and j = (-1) 1 / 2 . The other symbols are: r, 

radial distance; r 0 , the radius of the wire; and <5, the skin depth. In view 
of this distribution, calculations will be concerned mainly with wires whose 
diameter does not exceed 1 mm (number 18 AWG wire). 

The main parameter is the inductance which should be as low as possible in 
order to minimize the effects of switching. Consequently, emphasis will be 
placed on the determination of L which will be calculated two different ways: 
one method will utilize the definition of inductance (flux linkage per unit 
current); the other method will utilize the calculation of capacitance from 
which the inductance will be evaluated from the relation L = (Cv 2 ) -1 = ye/C 
= Z on /v6- 


2.2 Inductance Calculations 

In the following calculations, the permeability of all materials will be 
assumed to be that of free space (y 0 = 4irxl0~ 7 H/m) . The inductance per unit 
length of a bifilar lead of radius r a (current in each direction is assumed 
uniformly distributed over the cross section of the wire, see fig. 1) separated 
from the return lead of radius r a ' by a distance D aa ' is (ref. 8) 
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when the return wire radius is the same as in the first wire (the case here 
under consideration). For an n wire system (n pairs) one can form an n x n 
matrix 
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L , . . L L 

nl in nn 


where L-jj = Lj-j, and calculate the equivalent inductance of the system 
(refs. 9 and 10) 


( 3 ) 



H/m 


which can be determined to any given degree of accuracy (ref. 10). 


(4) 


The FORTRAN program for this calculation, CRIMPD.F, which is listed in 
appendix A, uses a series of menus to facilitate the entry of different two- 
dimensional conductor geometries. The entries are the coordinates of the 
centers of the conductors and their respective radii, and they may be entered 
either in rectangular or polar coordinates, (see figs. 2 and 3 for the polar 
coordinates) . 


2.3 Check on Calculations 

To check the results obtained from calculations performed according to 
the above discussion, two additional calculations were performed for several 
geometries. One check used the extension of equation (2) to multiple parallel 
conductors 
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with the parameters defined in figure 4. The equivalent inductance per unit 
length was calculated using 
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Calculations for 4-, 8-, and 16-conductor arrangements agreed exactly with the 
other calculations. 

The other check utilized the expression for capacitance between two 
infinitely long conductors of circular cross section (ref. 11) (radius r, 
distance between centers D) 



so that the external self-inductance of one conductor is 



because LC = ye . Calculations were performed using the above argument of the 
logarithmic function instead of using the D-jj's. The results, for the ratio 
of the distance between conductor centers to conductor radius equal to four, 
are about 7 percent lower than the other calculations. This was to be expected 
because the internal inductance, ti/8ir, does not appear in equation (8), but its 
contribution is negligible for the distance-to-radius ratio of four or more, 
and the 7 percent discrepancy is well within bounds (ref. 10). Note that the 
dimensional relations of this report meet this criterion: the ratio of the 

distance between centers to conductor radius never is less than four. Con- 
sequently, the error-if any-should not exceed 7 percent, a value usually 
acceptable in practice. 


3.0 CONFIGURATION DATA 
3.1 Circular Cylindrical Arrangements 

Calculations were performed for wire arrangements of the type shown in 
figures 2 and 3. For the single "ring" of wires, as in figure 2, it was found 
that the best arrangements (i.e., the lowest values of inductance and charac- 
teristic impedance) were for alternating current directions in adjacent wires. 
The values of inductance decrease with increasing number of wires. This cir- 
cular configuration was not suitable, however, because of the large number of 
wires needed to reduce the characteristic impedance to acceptable values (on 
the order of an ohm). Furthermore, with an increasing number of wires, the 
radius of the "ring" R is increasing, and this leads to inefficient use of 
space. 

A better configuration was of the multiring variety. Tables I and II list 
values of L per unit length, and of Z 0 for different number of wires (N 
represents the total number of wires in the line, with N/2 in the inner and N/2 
in the outer ring). The angle e is defined in figure 3 and represents the 
shift of the first wire in the outer ring relative to the first wire in the 
inner ring. As in the case of the single ring, the lowest values were obtained 
for alternating polarities. 
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One of the problems with the coaxial arrangement is that the inductances 
of the wires in the inner ring are lower so that more current is concentrated 
there relative to the outer ring. This can be compensated by reducing the 
radii of the wires in the inner loop. 

For comparison of the different wire radii, calculations for the arrange- 
ment shown in figure 5 were made. The values for the five ring arrangement for 
the same wire radii (Line 1), and for different radia (Line 2), are listed in 
Table III. These show that in the latter case, the overall inductance is 
increased which means that an increase in the number of wires would be required. 
Clearly, this again suggests not only an inefficient use of space, but also 
higher L and Z 0 values than comparable rectangular arrangement (see items 
1 and 2 of Table IV) . 


3.2 Rectangular Arrangement 

To carry 100 A, #4 AWG wire is recommended (ref. 12) whose cross section 
is 21.15xl0~ 6 m 2 . Inasmuch as the skin depth at 20 kHz is approximately 
0.5 mm, the largest reasonable wire would be #18 AWG (radius, 0.512 mm; cross- 
sectional area, 0.8231 mm 2 ). The minimum number of #18 AWG wires would there- 
fore be 26. In subsequent calculations, a larger number of wires will be used. 
(Note: the 26 wires do not include the return path which calls for additional 

26 wires. The minimum for the line would then be 52 wires #18 AWG.) 

Figures 6 to 8 show the results of 80-wire 4 by 20 rectangular arrangement 
which appears to have the best parameters: L = 10.4 nH/m; C = 2.5 nF/m; 

Z 0 = 2.1 ft (the computer printed values were, respectively, 10.4196, 2.5264, 
2.0611). As the method to compute these values was the "inductance" method, 
the numbers represent the upper bound. The wires in calculation were #18 AWG, 
total cross-sectional area in one direction is 32.68 mm 2 . 

A calculation for 256 wires #24 AWG in an 8 by 32 arrangement, as shown 
in figure 9, was also made. As expected, the inductance and characteristic 
impedance values are lower than for the 80-wire configuration: L = 3.2 nH/m, 

C = 7.9 nF/m, Z 0 = 0.6 ft (printout values, respectively, 3.2193, 7.9382, 
0.6368); cross-sectional area, 25.1 mm 2 . Note that the variations of the 
self-inductance of the individual wires are greater than for the 80-wire 
arrangement. 

Given the same size wires and the same distances among them, inductances 
of different size bundles are proportional to the inverse ratio of the number 
of conductors in the bundle. This was used to arrive at the values listed in 
Table V for the minimum number of wires required to pass 100-A currents, (if 
the self-inductances of each of the wires were the same, this relation would 
hold for any number of wires. Inasmuch as this is not the case (see figs. 7 
to 9), the number of wires should be greater than about 20 for the approxima- 
tion to be acceptable.) 

Table V lists several parameters for lines composed of a minimum number 
of conductors that can carry the rated (100-A) current. The line lengths are 
150 m each, so that the dc resistances were calculated for 300-m length to take 
into account the return path. The total number of wires in each line is twice 
that listed in the table. The corrections for the ac values were taken from 
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standard tables (ref. 13), and those for proximity effects from Smith 
(ref. 14). The proximity effects account for about additional 40 percent of 
resistance (see Table I and fig. 5 of ref. 14) for a distance-to-radius ratio 
of four (two in ref. 14 because the distance between the centers in this refer- 
ence is "2c"). 

To adjust the value of inductance per unit length and characteristic 
impedance of the 52-wire #18 AWG line listed in Table V, it is only necessary 
to multiply by 80/52 = 1.5, which would make 9.7 < L < 16 nH/m, and 1.9 < Z 0 
< 3.2 Q. Similar adjustments can be made for the other lines. 


3.3 Miscellaneous Considerations 

At the rated current, about 6.5 W will be dissipated per meter of line, 
which would mean an increase (ref. 15) between 0.08 to 0.13 K, depending upon 
the type of insulation. Such temperature rise is not sufficiency large to 
cause concern. 

The coefficient of expansion for copper is on the order of 10~ & m/K from 
25 to 1200 K (ref. 16). At 300 K, the expansion is 16.8xl0 -6 m, which would 
mean a total of 2.5 mm for the entire 150-m line. This, too, is not considered 
significant. 

For the recommended types of aircraft wire (ref. 17), Mil-W-16878, the 
dielectric strengths substantially exceed the specification of standing off 
1 kV, which - for the recommended configuration - would be on the order 
1 kV/mm = 1 MV/m = 25 V/mil. The breakdown strength of polyethylene exceeds 
20 kV/mm at 20 kHz, as does that of polystyrene. Teflon, and polypropylene 
(ref. 18). 

The conductance per unit length between the individual conductors can be 
calculated from the relation 


C/G = e/a (9) 

which, for the materials under consideration, will be on the order of 
10-15 m ho/m (R of the order of 1 0l 5 fi/m) , using 1 0^ 7 Q-m for resistivity 
(ref. 6). Clearly the losses due to transverse currents can be neglected. 


3.4 Connectors and Junction Boxes 

A conventional, circular, multipin connector is shown in figure 10 as an 
example. Junction boxes with connecting wires, with conventional connectors of 
circular or rectangular shape will fulfill the specifications. The junction 
box shown in figures 11(a) and (b) would be more suitable, however, because the 
interconnecting plates substituted for wires would reduce the inductance 
locally. Connectors could be of the conventional, multipin variety, with 
soldered or brazed connections internally. Alternatively, internally the pins 
could be connected to blades for contact with the conductors. The specific 
configuration is not very critical, but it would be recommended that the con- 
necting surfaces be gold-plated. The reason for this is to prevent increases 
of contact resistances as a result of fretting (ref. 19). 
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3.5 Miscellaneous 


In addition to the configurations described above, parallel-plate arrange- 
ments of the type shown in figures 12 and 13 were analyzed. That of figure 12, 
in spite of the very low characteristic impedance (0.14 Q) , is thought to be 
less suitable than the configurations recommended in the next section of the 
report because the thin plates will very likely buckle. Should the plates be 
made thicker, the line would be too heavy and too unwieldy. Also, connectors 
with suitable characteristics would be difficult to design. 

The program for the determination of line parameters for the configuration 
shown in figure 13, for different values of e, thickness, and spacing, is 
listed in appendix B. The values of inductances per unit length, and charac- 
teristic impedances are at least four times larger than those of the configura- 
tion shown in figure 6. 


4.0 RECOMMENDATIONS 

The line recommendations are made on the basis of the assumption that off- 
the-shelf availability, or ease for assembly-line production, are secondary to 
the other attributes such as low impedance, reliability, and the like. As a 
result, the best choices for the line are listed in Table IV. The total number 
of wires per 150-m line is listed, with one-half representing the return path. 
The direction of current alternates between neighboring lines. (The difference 
between Table IV and Table V is that in the latter one-half of the wires was 
used in a 300-m long loop for the determination of the total resistance, 
whereas the number of wires in the former represent two parallel sets each 
150 m long.) The alternating current directions have at least two advantages: 
lowering the inductance of the comparable wire arrangement with adjacent cur- 
rents in the same direction by at least a factor of four, and by reducing the 
external magnetic fields which reduces electromagnetic interference. Numbers 
1 and 6 in Table IV represent the configurations of figures 6 and 9; the 
remainder represents the minimum number of wires required for 100-A currents. 

Multipin connectors of the type generally available from manufacturers 
such as Litton, Amphenol, or others, would be satisfactory for the junction 
boxes shown in figure 11, with gold-plated contacts to prevent increases of 
resistance due to fretting (platinum plating would also be acceptable). 

A general recommendation for any type of line used in an aircraft power 
distribution system is that instead of one line 150 m long there should be 
several 5-, 10-, or 15-m sections connected in series via junction boxes. The 
advantage of this would be that in case of failure of a section, it could be 
easily replaced-even in flight. Also, the several junction boxes at given 
intervals would facilitate connections to the respective loads. 
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PROGRAM C R IMP D.F 


DOUBLEPRECISION A,E,C 

DOUBLEPRECISION SCI : 1024) ,0(1 : 1 024),X (1 : 1 024) , Y (1 :'1024) 
DOUBLEPRECISION L (0 : 1 024) ,L0 CO: 1 024 ) , LI (0: 1024),LN 
DOUBLEPRECISION D (0: 1 024) , T (1 : 1 024) , R ( 1 : 1 024) , B <0: 1 024) 
DOUBLEPRECISION S A/ S D, P I, S LOPE ANG 
INTEGER M,K,UUU,JE,ALTPOL 
INTEGER U,J,I,N,RP(1:1Q24) 

CHARACTER Q1 , POL( 1 : 1 024 ) * 3 
PARAMETER (P 1=3 .1 41 592654) 

Q1 =* SELECTION IS OUT OF BOUNDS. TRY AGAIN...' 

CONTINUE 
DO 140 1=1,64 
PRINT* 

S(I)=0. 

0(I>=0. 

X(I)=0. 

Y(I ) =0. 

L ( I ) = 0 . 

Loa)=o. 

LI (I)=0. 

o(i)=o. 

TCI)=0. 

R 1 1 ) = 0. 

B ( I ) =6. 

POLE I ) = ' N ’ 

CONTINUE 

N=1 

PRINT’ (////////)' 

PRINT*," Main Menu' 


PRINT* 

PRINT*,' 

PRINT*, ' 

PRINT*,' 

PRINT*,' 

PRINT*, * 

PRINT*,' 

PRINT*,' 

PRINT*,' 

PRINT' ( 7 / ) ' 


1. Polar coordinate input* 

2. Rectangular coordinate input' 

3. Data file input' 

4. Run inductance calculation (Skilling method)' 

5. Run capacitance calculation (Cheng method)' 

6. Check individual conductors' 

7. Exit to shell* 

3. Enter new configuration (reset index counter) 


PRINT*, "Choose from 1 thru 8...' 

R E AD*,U 

I F (U .EQ. 7) GOTO 7000 
I F (U .EQ. 6) GOTO 6000 
I F (U .EQ. 5 ) GOTO 250 
I F (U .EQ. 4) GOTO 4000 
IF (U .EQ. 3) GOTO 250 
I F (U .EQ. 2) GOTO 2000 
I F (U .EQ. 1 ) GOTO 1000 
I F (U .EQ. 8) GOTO 100 
PRINT*, Q1 
GOTO 200 

PRINT*, 'ROUTINE NOT OPERATIVE, CHOOSE ANOTHER MODE... 
GOTO 200 


Polar coordinate input 

PRINT*,' Polar coordinate input' 

PRINT*,' 1. Automatic ring program' 

PRINT*,' 2. Individual input’ 

PRINT*,' 3. Exit to main menu' 

PRINT*, 'Choose one...' 

READ*,U 

I F (U .EQ. 3) GOTO 200 
I F (U .EQ. 2) GOTO 1500 
I F ( U .EQ. 1)G0T0 1100 
PRINT*, Q1 
GOTO 1000 

Auto polar ring 

PRINT*, 'AUTO POLAR RING’ , . , 

PRINT*, 'Input number of conductors in this ring...' 
READ*,M 

SA=2.0*PI/REAL(M> 

SD=SA*1 80.0/PI 



APPENDIX A. - PROGRAM FOR INDUCTANCE CALCULATIONS 



* * * * 


1110 


1120 


1160 


1180 

1200 

1220 

1260 


1230 

1300 


1320 

1340 


1500 


PR INT*/ ' Input coordinates/ radius and polarity of 1st conductor* 
PRINT*/*in this ring Ir (mm) / 1 he t a (degr ees ) , raai us turn) / + or ->• 

RE AD (5/ 11 10) SCN)/TCN)/RCN)/PQLCN) 


F0RMATCF8.4/F8.4/F8.4/A1) 

OCN)=TCN)*PI/l80.0 

SCN)=SCN)/1000.0 

RCN)=RCN>/1000.0 

X(N)=SCN)*C0S10CN)) 

YCN)=SCN)*SINCOCN)> 

RP <N) = 0 


PRINT* (///A)'/* Choose polarity of wires... 1 
PRINT*/' 1. Alternating* 


2 . 

1: 


All the same as 1st' 
Individually selected* 
Exit to last menu* 


11 GOTO 
21G0T0 

3) GOTO 

4) GOTO 


1200 

1220 

1180 

1000 


1560 


PRINT*/' 

PRINT*/ * 

PRINT*/ * 

RE AD*/U 
I F (U -EQ. 

IFCU .EQ. 

IFIU .EQ. 

IFCU .EQ. 

PRINT*, Q1 
GOTO 1120 

UUU=1 
GOTO 1260 
UU=0 

GOTO 1260 
UU=1 

DO 1340 I=N+1,N+M-1 
RP(I)=0 
S(I)“S(N) 

0(11=0(1-1) + SA 
TCI)=T(I-1) + SD 
XCI>=S<I>*COSCOCI>) 

YCI)=SCI)*SINCOCI)) 

R(I)=RCN) 

IFCUUU .NE. 1 ) GOTO 1300 

PRINT* </,T2,A,I3,A,F8.5,2X,F8.4) • / 1 C onduct or * / 1 /'with coor ds . * , S ( I) *1 000. 0,T (I ) 
PRINT*/ ' Input polarity ( + or -)...' 

READ' (A) * / POLC I ) 

IFCU .EQ. 21G0T0 1320 
IFCP0LCI-11 .EQ. *+* )P0LCI)=’-' 

IFCPOLCI-1) .EQ. •-' ) POLC 11 = * + * 

GOTO 1340 
P0LCI)=P0LCI-1) 

CONTINUE 

N=N+M 
GOTO 200 

Individual polar input 

PRINT*/' Individual Polar Inputs* 

PRINT*/*Hov many conductors in this series?* 

RE AD*/K 

DO 1560 I=N,N+K-1 

PRINT*/ ' Input r Imm), t he ta Cde gre es ) , rad i us of conduct or C mm)/ po lari ty C+ or -)...* 
READ (5/1110) SCI)/TCI)/RCI)/POLCI) 

RP C I ) =0 

RCI)=RCI)/1 000.0 
SCI)=SCI)/1000.0 
0CI)=TCI)*PI/180. 0 
XCI)=SCI)*COSCOCI)) 

YCI)=SCI)*SINCOCI)) 

CONTINUE 
N=N*K 
GOTO 200 


2000 

2100 


rectangular coordinate input 
CONTINUE 

PRINT*, 'RECTANGULAR COORDINATE INPUT* 

PRINT*/ * Choose one...' 

PRINT* 

PRINT*,' 1. Individual input' 

PRINT*/' 2. Line input Cequally spaced wires) 1 

PRINT*/* 3. Return to main menu* 

READ*,U 

IFCU -EQ. 1 ) GOTO 2120 
IFCU .EQ. 2) GOTO 2200 
IFCU .EQ. 3) GOTO 200 



2120 


2H0 

2150 


* 

* AUTO 
2200 


221 0 
2230 


2300 

2320 


PRINT*/' Individual Rectangular Inputs' 
PRINT*/'How many conductors In this series-..' 

RE AO*/ K 

DO 2150 I-N/N+K“1 
PRINT* 

PRINT*/ ' Input coordi n at e s/radi us & polarity of 
PRINT*/ ' (X C mm) / Y(ram) / r< mm) /+ or -)* 

RE AD (5/1 110) X(I)/Y(I)/RlI)/P0L(I) 

RP(I)=1 

X(I)=X<I)/1000.0 
Y(I)=Y(I)/1 000.0 
R(I)=R(I)/1000.0 
S(I)=SQRT (X (I)**2«0*Y (I) **2.0) 

IFCS(I) .EQ. 0.0)G0TO 2140 
T(I)=180.0/PI*AC0S{X(I)/S(I)) 

GOTO 2150 
T ( I ) =0 . 0 
CONTINUE 
N=N+K 
GOTO 200 

LINE 

CONTINUE 

PRINT*/'AUTO LINE INPUT* 

PRINT*/ ' Input coordinates/ radius & polarity of 
PRINT*/'CX(mm)/Y(mm)/r(mm)/+ or -)...' 

READ (5/1 110) X(N),Y(N)/R(N)/POL(N) 

X(N)=X(N)/1 000.0 
Y(N)=Y(N)/1 000.0 
R(N)=R(N)/1 000.0 

S(N)=SQRT(X<N)**2.0 + Y(N)**2.0> 

I F ( S ( N) .EQ. 0.0) GOTO 2210 
T(N)=180.0/PI*AC0S(X(N)/S(N)) 

GOTO 2230 
T (N) =0.0 
CONTINUE 
RP(N)=1 
PRINT* 

PRINT*/ ’ Input number of conductors on this line. 
READ*/M 

PRINT*/ * Input distance between centers of conds 
READ*/SPACING 

SP AC I NG= SPACING/1 000.0 

PRINT*/ ' Input slope angle of the line (positive 
READ*,SLOPEANG 

SLOPE ANG=S LOP EANG*P 1/1 80.0 
DO 2300 I=NFl/N+M-1 

X(I)=X(I-1) + SPACING*COS(SLOPEANG) 

Y(I)=Y(I-1) + SPACING*SIN(SLOPEANG) 

R(I)=R(I-1) 

S(I)=S9RT(X(I)**2+YCI)**2) 
T(I)=ACOS(X(I)/S(I))*18 0.0/PI 
RP ( I ) =1 
CONTINUE 
CONTINUE 

PRINT' (///A) '/' Choose polarity of wires...' 

" ‘ 1 . Al t erna t ing ' 

2. All the same as the 1st* 

3. Individually selected* 

4. Return to last menu' 


2340 


PRINT* 
PRINT*. 
PRINT*/' 
PRINT*/' 
READ*/U 
I F ( U .EQ 


.. . _ . 1 )GOTO 2340 
I F ( 0 .EQ. 2 ) GOTO 2380 
IF (U .EQ. 3 ) GOT 0 2420 
EQ. 4 ) GOT 0 2000 


I F ( 0 

PRI NT*/ Q1 
GOTO 2320 
CONTINUE 

IF(P0L(N) .EQ. ' + ' ) T H EN 
ALTP0L=1 

ELSEIF(POL(N) .EQ. ' - * ) THEN 
ALTPOL=-1 
ENDIF 

DO 2360 I=N+1/N+H-1 

I F ( ALTPOL .EQ. 1)THEN 
POL ( I) = '”* 

ELSEIF (ALTPOL .EQ. -1 ) THEN 
P0L(I>='*' 

ENDIF 

ALTP0L=“ ALTPOL 


cond. ' / 1/ ' . * 


1st conductor...' 


on this line (SPACING) . 
x-axis=0 degrees)' 



2360 

CONTINUE 
N=N*M 
GOTO 200 

2380 

CONTINUE 

00 2400 I=N+1,N*M-1 
POL (I ) =POL( N) 

2400 

CONTINUE 
N=N+M 
GOTO 200 

2420 

CONTINUE 

DO 2440 I=Nt1/N*N-1 


PRINT" (//T2/A/I3/A/F8.5/2X/F8. 4)' /"conductor"/!/' with coords. 
PRINT*/ ' Input polarity ( + or 
READ*/ POL (I ) 

2440 CONTINUE 

N=N+M 
GOTO 200 


x(mm)="/X(I)/"y(mm)="/Y(I> 


3000 


4000 


4050 

4140 


4160 

4180 

4200 


4220 


4240 

4260 

4280 

4300 

4320 


4400 

4420 

4440 


4460 


cont i nue 


D ( 0) =1 . 0 

PRINT*/ * Input value of dielectric constant..." 

R E A0 * / E 

print* 

00 4320 I=1/N-1 

IF (POL (I ) . EQ . 1 - • ) GOTO 4220 
00 4200 J-1 /N“1 

IFU .EQ. I > GOTO 4160 

DCJ)=SQRT((X(J)-X(I)>**2+(Y(J)-Y(I))**2) 

IF (P0L(J ) .EQ. "-* ) GOTO 4180 
0(J)=1 .0/0(0) 

GOTO 4180 
D(J)=1 .0 

D(J)=D(J)*0(J-1) 

CONTINUE 
GOTO 4300 

0 ( 0) =1 . 0 
00 4280 J=1/N-1 

I F ( J .EQ. I ) GOTO 4240 

0CJ)=SQRT((X(I>-X(J))**2+(Y(I)-Y(J))**2) 
IF(POLIJ) .EQ. • + • ) GOTO 4260 
D(J)=1 -0/D(J) 

GOTO 4260 

o(j)-i .o 

D(J)=D(J)*D(J-1) 

CONTINUE 

L(I)=(0.5+2.0*L0G(0(N-1 )/R(I)) ) *100.0 
CONTINUE 


LO (0) =0. 0 
LI (0)=0.0 
00 4440 I=1/N-1 
IF (POL(I) .EQ. 
LI (I)=1 .0/L(I) 
LO ( I) =0.0 


' - • ) GOTO 4400 


GOTO 4420 

‘ 1 = 0.0 


LI (I)= 

L0(I)=1.0/L(I) 

LI ( I ) =L1 (I)+L1 (1-1 ) 

L0(I)=L0(I)+L0(I-1) 

CONTINUE 

B(0)=0.0 
00 4460 1=1 / N~1 
8(I)=R(I)*1 000.0 
A=(B(I))**2+(B(I-1))**2 
3(I)=SQRT(A) 

CONTINUE 

A=A*PI/2.0 

I F ( ASS ( LO ( N~1 ) ) .L£. 0. 0000001 ) GOTO 4480 
LN=1.0/L1 (N-1 )+1.0/LO(N-1 ) 



* * 


4470 goto 4490 
4480 LN=1.0/L1(N-1) 

4490 CONTINUE 

C“E *100.0/9. 0/LN 

PRINT' (T2/A/F10. 4, A)', 'INDUCTANCE OF CIRCUIT IS L= '/LN/ 'nH/m' 

PRINT'(T2/A/Fj0.4/A)*/'CAPACITANCE IS C='/C/'nF/m' 

PRINT* (t2/a/F10.4/a)'/’CHAR. IMPED. IS I=* / SORT (LN/C) / ' ohms ' 

PRINT* (T2/A/T51/F10.4/T62/A//) '/'TOTAL CROSS-SECT. AREA IN ONE DIRECTION IS A='/A/'sq. ram' 

PRINT'!//)' 

t 

PRINT’ (2X/A/T10/A/T20/ A/T 30/ A/T 37/ A/T50/ A/T 60/ A/T72/ A) ' / ' eond' / • x (mm) ' / * y ( ram) ’/'r (ram) '/■'theta (deg) '/'polarity'/ 'radius '/■'l' 

00 4640 1 = 1 / N-1 

IF(RP(I) .EQ. DGOTO 4600 

PRINT' (2X/ I3/T2 5/ F10. 4/ T36/ FI 0. 4/TS3/ A/T56/ FI 0.4/T67/ FI 0.4) ' / 1/ S (I) * 1 000.0/ T( I) /POL (I )/R (I )*1 00Q.0/LII) 

GOTO 4640 

* 

4600 continue 

PRINT' (2X/I3/T7/F10.4/T18/F10.4/T53/A/T56/F10.4/T67/F10.4)'/I/X(I)*1000.0/Y(I)*1000.0/POL(I)/R{I)*1000.0/L(I) 

4640 CONTINUE 
4700 PRINT* 

goto 200 

It 

lr 

5000 continue 


6000 


6200 

6300 


PRINT*/' Check individual conductors' 

PRINT*/ ' Input number of 1st conductor to be checked...' 

RE AD*/ J 

PR INT*/ * Input number of last conductor to be checked...' 

RE AO * / J E 

PRINT' (2X/A/T9/A/T18/A/T33/A/T43/A/T50/A/T60/A) '/' cond' /'r(mm) •/• theta(deg) '/'x(mm) ' /'y( mm) '/'polarity'/' radius* 

DO 6200 I=J/JE 

PRINT* (2X/I3/T7/F8.4/T21/F9. 3/ T31 / F3 . 4/ T41 / F 8. 4/T5 4/A/T 58/ F8. 4 ) * / 1/ S ( I ) *1 000. 0/ T( I ) / X ( I ) *1 000. 0/ Y ( 1 ) *1 000 . 0/ POL (I ) / R (I ) * 1 000 
CONTINUE 

PRINT*/ ' Check another conductor? (1 =y e s/0=no) . . . * 

RE AD*/ U 

IF (U .EQ. DGOTO 6000 
I F (U .EQ. 0) GOTO 200 
PR INT*/ Q1 
GOTO 6300 


7000 


7500 


PRINT*/'Are you sure you want to exit proqram?’ 
PRINT*/' 1 — YE S' 


PRINT*/ ' 
RE AD*/U 
I F (U .EQ. 
IF (U .EQ. 
PRINT */Q1 
GOTO 7000 
PRINT* 

END 


2=N0 ' 

2) GOTO 200 
DGOTO 7500 



PROGRAM 
*speci f i cati ons 
INTEGER 
INTEGER 
INTEGER 
INTEGER 
INTEGER 
REAL 
REAL 


LAPLAC 


NMIN, NMAX, MM IN, MM AX 
P(0: 64,0: 64),L- 
N,M,N1,N2,N3,N4,M1,M2,NS,MS 
LMAX , LS,OUT,IN,ID,JD,LT 

I,J,K,K0,KQ,KT,KMAX,JSL,MM2,MM3,M2SL,RPN1,RPM1 
SL,ERLIM,ACAP,DCAP,RPA,RP9,RPC 
ERMAX, LAMBDA 


REAL F(0:64, 0:64), C(G:18, 1:5), H,RP1 ,RP2 
CHARACTER PS(0:18)*1 
COMMON F,P,C,L,H,RP1,RP2 
COMMON/ REL AX /LAM8DA,ERMAX,I,J 
CO MMON/CONTQR/NMI N, NMAX, MMIN, MMAX 

* i n i t i alizations- — - — 

DATA N, Mr NMI N,NMAX, MMIN ,MMAX /3 2, 64, Or 20*0*43/ 

DATA ERLIM,KMAX/1E~4,100/ 

DATA N1,N2,N3,N4,M1,M2,N$,MS/4,12,16,4,16,32,6,18/ 
DATA (PS(I),I=0,1 8)/ 


'•'•b'/'e'/'fl'/'eVf • 


DATA LAMBDA,LMAX,H/1.4,4,0.1 25/ 

DATA RP1 ,RP2, RPN1 ,RPM1 , XD, KT / 1 ,1 , 1 6, 32, 0, 0/ 

♦ function statements 

I RND ( FI >=NINT(F1*1E3) 

♦ data Inputs—-—-- — 

PRINT* (A) ', 'OUTPUT TYPE <0=ALL, 1 = FINAL, 2=DAT A FILE)?’ 
READ*, OUT 

PRINT* (A)*, * INPUT TYPE (0 = INTERNAL, 1 =EXTERNAL, 2 = D ATA FILE)* 
READ*, IN 

IF (IN .EG. 2) THEN 

READ* <3I3,2F6.3)*,N,M,L,H,SL 
DO 2 JD=0,M 
DO 1 I D=0,N 

READ* (2 13, El 5. 7, 1 3) ' , I, J , F ( ID, J D) ,P C I D, J D) 

1 CONTINUE 

2 CONTINUE 
CALL CNTOUR 
GOTO 999 

ENDIF 

I F (IN .LT. 1)THEN 

PRINT* (A)', 'INPUT SLOPE ( 0,0. 5, 1 , 2,3=NO FLAP ):* 

READ*, SL 
ENDIF 

3 PRINT* (/T2,A,F6.3,A//,10A5,2A8)*, 

: 'BOUNDARY DI MEN SI ON S ( 1 = * ,H , * mm ):*, 

: 'N','M', , N1*,'N2*,*N3','N4*,'M1*,*M2 , ,*NS','MS', , SL0PE*, , LA 

IF (IN .LT. 1 ) GOTO 4 

READ*,N,M,N1,N2,N3,N4,M1,M2,NS,MS,SL,LAMBDA 

4 PRINT' (10I5,2F8.3//)',N,M,N1,N2,N3,N4,M1,M2,NS,MS,SL,LAMBDA 
PRINT'(T3,A//,4A6)', 

: 'INPUT CONTOUR LI MITS* , * NMIN * , ' NMAX * , ’ MMIN * , * MMAX * 

IF (IN .LT. 1 ) GOTO 5 
RE AD*, NMIN, NMAX, MMIN, MM AX 

5 PRINT* (416/) NMIN, NMAX, MMIN, MM AX 


h*,*i*,'j*,* 


l*/ 


APPENDIX B. - PROGRAM FOR "STRIPLINE 11 CALCULATIONS 



♦boundary conditions- 
PRINT* (A)* ‘ 


0)THEN 


1 r * OUTSIDE LAYER THICKNES$=?' 


. HOMOGENEOUS 

READ^RPl 
IF (RPl -GT. 

KO = 0 
RP2=RP1 
RPl =1-0 
ELSE 

RP2=A3S(RP1 ) 

RPl-1.0 
K0 = 0 

PRINT* (A) * - * DIELECTRIC 
READ** KD 

IFCKD -EQ- 1)THEN 
KT=0 

PRINT* (A)’#' DISTANCE FROM C ENT ER = ? * 
READ**KT 
RPM1 =M1 +KT 
RPN1 =N3 
ELSE 

PRINT' (A) * 

R£AD**KT 
RPN1 =N3+KT 
RPM1 =M2+KT 
END I F 
ENDIF 
RP A=2*RP1 
RP8=2 *RP2 
RPC=RPl +RP2 

♦point type array setup 

♦homogeneous Interior point 

DO 7 J=1/5 
C (5* J ) = RPA 

7 CONTINUE 
DO 8 J = 1 r 5 

C<6* J)=RPB 

8 CONTINUE 

♦dielctrlc In two adjacent quadrants 
C(7*1)=RPC 
C(7*2)=RPC 
C(7*3)=RPA 
C(7*4)=RP8 
C(7*5)=RPC*2 
CC8*1 )=RPA 
C(8*2)=RPB 
C ( 8*3) =RPC 
C<8#4)=RPC 
C ( 8*5 ) = RPC *2 
C ( 9*1 ) = RPC 
C<9*2)=RPC 
C(9*3)=RPB 
C(9*4)=RPA 
C(9* 5) = 2*RPC 
C<10*1)=RPA 
C < 1 0*2) =RPB 
C ( 1 0*3) =RPC 
C<10*4)=RPC 


DIELCTRIC C0NSTANT(YES=GE(1)/N0=LE(-1)> 


CONFIGURATION ( 1=IN.$IDE* 2*QUT SID E ) * 



. C <1 Or 5) =2*RPC 

♦dielctric In one quadrant - "— — 
C(11r1)=RPA 
C(11r2)=RPC 
C ( 1 1 r 3) =RP A 
C£11r4)=RPC 
C(11r5)=RPC+RPA 
C ( 1 2# 1 ) -RPA 
C(12r2)=RPC 
C(12#3)=RPC 
C(12r4)=RPA 
C112,5)=RPC+RPA 
C ( 1 3r 1 ) =RP A 
C ( 13r 2) =RPC 
C(13/3)=RPA 
C ( 1 3r 4) -RPC 
CC13r5)=RPC+RPA 
C ( 1 4r 1 ) -RPC 
C(14r2)=RPA 
C£14r3)=RPA 
C £ 1 4r 4) -RPC 
, C C 14r 5) =RPC*RPA 

♦ dielectric In three quradrants- 

C(15r13=RPC 

C<15,2)=RPB 

C<15r3)=RPC 

C Cl 5r4 ) -RPB 

C(15r5)=RPC+RPB 

C C 1 6r 1 ) =RP C 

C £ 1 6r 2) =RPB 

C(16r3)=RP8 

C(16r4)=RPC 

C(16r5)=RPC+RPB 

C£17r1)=RPC 

C£17r2)=RPB 

C£17r3)=RPC 

C£17r4)=RPB 

C(17r5)=RPC+RPB 

C(18r1)=RP8 

C < 1 8r 2) = RPC 

C(18r3)=RPC 

C (1 8r4) =RPB 

CC18r5)=RPC+RPB 

♦ Initial point type — 

DO 10 I “ 0 r N 
DO 9 J=0 ,n 
F£IrJ)*0,0 
DC ( Ir J )-RPl 

9 P ( I r J ) *5 

10 CONTINUE 

♦defined boundaries 

DO 600 LS=2r0r-1 
L=£2**LS) 
KQ=JCMAX/R£ALIL) 

♦upper boundary (xry-m)* * 

DO 15 I=0#NrL 
P£IrM>=2 




15 

*ri ght 
20 

* lower 


30 


35 

♦ left 


40 

♦def In 
♦cente 


50 

55 

♦outer 


60 

65 

♦outer 


70 

75 

♦secon 

80 


120 

140 


CONTINUE 

boundary (x=n/y)* 

DO 20 J=0*M,L 
P(N,J)=2 
CONTINUE 

dielectric boundary (x*y=0)* — 

DO 30 I = (N1+L)/- (N2-L)xL 
P(I,0)=3 
CONTINUE 

DO 35 I=(N3+L),N,L 
P ( I / 0) =3 
CONTINUE 

dielectric boundary (x=0*y)* 

DO 40 J=(Ml+L)rNxL 
P(0 aJ)=4 
CONTINUE 

ed conductors — 

r conductor (x*y)* — 

DO 55 I=0,N1,L 
DO 50 J=0,M1*L 
F(I#J)*1.0 
P<I,J)=1 
CONTINUE 

conductor (KrV ) * 

DO 65 I=N2,N3,L 
DO 60 J—Or M2^L 
P(I/J)=2 
CONTINUE 

conductor edg e ( x ✓ y ) * --- 

IF(SL -GT. 2) GOTO 80 
HM2=M2-LMAX 
MM3=H2+5*LHAX 
DO 75 J=MM2,«M3,L 
DO 70 I=N3* N4j--L 

JSL ss NINTlSL*(N2~I)/L)*L 

M2$L=M2+JSL 

IF(((M2SL-<N3-N2+L>) .LT. J -AND, J -LE. M2SL) 
(N2 -LE. I .AND, J .LT. M2SL) ) THEN 
F (I* J)*0.0 
P (I* J ) =2 
ELSE 

P(Ir J ) =5 
ENDIF 
CONTINUE 
CONTINUE 
d dielectric 
IFCKD .GT. 


fill 

. OJTHEN 

DO 140 d=1,CRPM1-L) 

DO 120 I*1r(RPNl-L) 

IF(P(I#J) .GT. 4) THEN 
P(I, J)=6 
ENDIF 
CONTINUE 

IFCKD .EQ. 2)THEN 
P CRPN1 * J ) =8 
ENDIF 
CONTINUE 


■ OR. 



160 

*air poc 
170 


★point t 
180 


190 

200 

★lambda 

★iterate 


4)THEN 


DO 160 I=1,(RPN1~1) 
IF(P(I,RPM1 ) .GT. 

PCI,RPM1>=7 
END IF 
CONTINUE 

IF (KD .EQ. 2) THEN 
P(RPN1,RPM1)=11 
ENDIF 
6NDIF 

ket Location 

IF (L .GT. 1 ) GOTO 180 
PRINT* ( A) * INPUT AIR 

READ*, I, J 
LT = 1 

I F ( I .GT. 0) THEN 
P(I,J)=5 
P < I~LT, J ) = 5 
PCI+LT, J )»5 
P(I,J-LT)=5 
P(I,J+LT)=5 
PCI-2*LT,J)=7 
PU,J-2*LT>=8 
P(I+2*LT,J)=9 
P(I,Jt2*LT)*10 
P ( I“LT, J**L) =1 1 
PCI+LT , J ~L) =1 2 
P ( I+LT, J+L) =1 3 
PC I~LT , J+L) =1 4 
P (I“2*LT ~L) =1 
PCI-LT,J-2*L)=1_ 

PCI+LT, J-2*L)=16 
P(I+2*LT,J-L)=16 
P(I+2*LT,J + L)-17 
PCI+LT,J+2*L>=17 
P(I~LT,J+2*L)=1 8 
P(I~2*LT,J+L)-1 8 
GOTO 170 
ENDIF 

ype printout — 

IF (OUT .EQ. 1 .AND. L 
IF (OUT .GT. UGOTO 200 
PRINT* v/T2,A//,T5,A,I2,/>*, 
* BOUNDARY CONDITIONS*^ * MESH 


POCKET LOCATION!! < I, J ) , <-1 ,X) = NONE 3 * 


5 

5 


GT. 1 ) GOTO 200 


SIZE = ' ,L 


PRINT* CT5,A/,T5,A/)*,*+:V=1000, *:V=0* ,* J : dV/dy=0, -:dV/dx~0* 

DO 190 J~M,Q,“L 

PRINT* <T6,33A)*,<PS{P<I,J>),I=0,N,L> 

CONTINUE 
PRINT* (////) * 

CONTINUE 

loop 

solution * 

DO 230 K=1,KQ 
IFCK . EQ* 1) 

ERMAX-1 
ELSE 
£RMAX=0 
ENDIF 


THEN 



210 

220 

230 

240 

250 


♦output 

300 


240 


320 


00 220 J — 0/ M/ L 
00 210 1 = 0/ Hr L 
CALL FT CALC 
CONTINUE 
CONTINUE 

IF( ERMAX .LE. ERLIM) GOTO 
CONTINUE 
K“0 

CONTINUE 

IF (L .LE. 1 .OR. OUT .LT. 1 ) THEN 
CALL CAPCT(NS/M$/ACAP/DCAP) 

END IF 

IfIoSt^EQ. 2)G0T0 500 

IF (L -GT. 1 .AND. OUT . GT. 0)G0T0 

IF (K .EQ. 0) THEN 

PRINT* (T2/A/I4/A //)*/ 

'SOLUTION HAS NOT CONVERGED IN 
ELSE 

PRINT* (T2/A/I4/A //>*/ 

•SOLUTION CONVERGED AFTER 
ENDIF 

PRINT* (T2/A/F9.6/Z) ' 

•WITH 


600 


'sKQj 


ITERATIONS 




*/K/* 


ITERATIONS* 


: / f\r I 

MAXIMUM 


»/ / t / 

ABSOLUTE 


ERROR OF* /ERMAX 


PRINT* (T2/A/I3/A///T2/A/F6.3///) '/ 


MESH SIZE=*/L/* 

•i * 


00 


PRINT* (T2/ A/ F1Q.6/A/)*/* CORNER 
PRINT *(/////////////)» 

♦ solution array printout 

IFCLS «^.T. _1 ) THEN 


♦ STEP SIZE*/ * LAMBDA = ' /LAMBDA 

CAPACITANCE (AIR) * * / AC A P/ * p F/m ' 


400 


♦data f 
500 


PLOT* 


/* t : v=1000 


* : V=0 


: V<1 O' 


510 

520 


♦terml n 
600 


PRINT* (T2/A/) '/'SOLUTION AT n*m POINTS/ <x=horiz/ y=vert) 
DO 400 J=M/0/-L 

PRINT' (T2/33(I4))*/(IRND(F(I/J))/I=0/N/L) 

CONTINUE 

PRINT 4 (////T2/A/7/T2/A//) * 

'CONSTANT VOLTAGE CONTOUR 
CALL CNTOUR 
ENDIF 
GOTO 600 

le output 

I F CL .EQ. 1) THEN 

PRINT* (3I3/2F6.3) '/N/M/L/H/SL 
DO 520 J=0/M 
DO 510 I=0/N 

PRINT* (213/ El 5.7/13) ' /I/ J/FCI/J)/ PCI/ J) 

CONTINUE 
CONTINUE 
PRINT* CI3) »/-1 
GOTO 999 
ENDIF 

ati on 

CONTINUE 

PRINT* (//A/zAz////) ' /' 

* 2=LAM8DA/ 3=B0UNDARY 
READ* (13) */I 
I F (I .LT. 01G0T0 999 
GO TO (62 0/630/ 640) /I 


. • 


CHANGE PARAMETERS? 
CONDITIONS) * 


(“1 TO END/ 1 = FLAP SLOPE/*/ 



620 RE AD*x SL 

GOTO 3 

630 READ**LAMBDA 

GOTO 3 
640 IN-1 

GOTO 3 
999 END 

************************************************************************* 
SUBROUTINE PTCALC 

♦ specif 1 cations — 

INTEGER I* J*PK 
INTEGER P(0:64*0:64)*L 
REAL F1*F2*F3*F4*F5*F6*F7*F8 
REAL C1*C2*C3*C4*C5 

REAL PAST*TEMP*ER*FINT* FNORM* FOI EL 
REAL ERMAXxLAMBDA 

REAL F(0:64*0:64)*C(0:18*1 :5)*H*RP1*RP2 
COMMON FxPxCxLxHxRPl*RP2 
COMMON/RELAX/LAMBDAxERMAXxI* J 

♦ function statements 

FI NT (FI *F2*F3*F4*F5*F6*F7*F8)=(4*(F1+F2+F3+F4)+(F5+F6+F7+F8))*0.05 
FN0RM(F1*F2)=(4^F2-F1)/3 

FDIEL(F1*F2xF3xF4*C1xC2xC3xC4xC5)=(Cl*F1+C2*F2+C3*F3+C4*F4)7(2*C5) 

♦ relaxation * — 

PK-P ( I* J ) 

GOTO (195*195*120, 1 40*1 80* 1 80* 1 60* 1 60* 1 60* 1 60) *PK 
120 TEMP=FN0RM(F(I*(2*L))xF(IxL)) 

GOTO 190 

140 TEMP~FNORM( F ((2*L)*J)*F(L*J)) 

GOTO 190 

160 TEMP-FDI EL (FC(I+L)*J)*F((I-L)*J)*F(I*(J+L))*F(I*(J“L))* 

: C(PK*1)*C(PK*2)*C(PK*3>*C(PK*4)*C(PK*5>) 

GOTO 190 

180 TEMP=FINT(F(I+L*J)*F(I~L*J)*F(I*J+L)*F(I* J“L>* 

^ : F(I+L*J+L)*F(I*L*J-L)*F(I-L*J+L)*F(I-L*J-L>) 

190 PAST=F(I*J) 

F ( I* J)=PAST+LAMBD A* (TEMP-PAST) 

ER-A8S(PAST“F (IxJ)) 

ERMAX=MAX(ERMAX*£R) 

195 RETURN 

END 

A ★ ★ ★ ★ * ★ * * it * *■ * * A it it * * k * ik -k "k 1c i* * * ★ ★ * ★ * it k Ir ★ * * * * * it ★ if he ★ * * k if -k A ft * k jfr * ★ j*r ★ Jr * A It it k ik irkic A 


SUBROUTINE CAPCT(NS*MS*ACAP*DCAP) 

♦ specificati ons 

INTEGER NS* MS 

INTEGER P (0: 64*0: 64) *L 

REAL ENORM*PERMFS*ACAP*DCAP 

REAL F( 0:64* 0:64) *C (0:1 8* 1 : 5 ) *H*RPl *RP2 

COMMON F*P*C*L*H*RP1*RP2 

PARAMETER (P£RMFS=8. 85419) 

♦ surface integration — 

ENORM-O.O 

ENORM=(F(NS+L*0)-F(N2-L*0)+F(NS+L*MS)-F(NS-L*MS))*0.5 
EN0RM=EN0RM>(F(0*HS+L)-F(0*MS-L)+F(NS*MS+L) -F(N5*MS-L>) *0.5 
DO 10 J*L*(MS-L)*L 

ENORM=ENORM+F(NS+L*J)-F(NS“L*J) 

10 CONTINUE 



, 0 l # , 1'/'2 , / , 3V4 , f , 5 , / , 6'/ , 7 , /'8 , / , 9'/ , 0‘/ 


DO 20 I~L* (NS-D/L 

EN0RM=EN0RM+F<I,MS+L)-f <IrMS~L> 

20 CONTINUE 

ACAP=A8$((5E-4*PERMFS*(NS+«S)*EN0RM)/H) 

DCAP=RP2*ACAP 

RETURN 

END 

SUBROUTINE CNTOUR 

♦ sped f 1 ca ti ons 

INTEGER IfirIC>-NCOL*MROW*MARKX/’MARKY 
INTEGER l,J,HQ,HQ,OH,DH,PO 
INTEGER NMINrNMAX,MMINrMMAX 
INTEGER P(0:64,0:64),L 

REAL F<0:64r0:64),C<0:18,1:5),rH,RPl,RP2 

CHARACTER D AR Y (0: 13 2,0: 260) *1 ✓ C S < 0: 1 6) * 1 

COMMON f,P,C,L,H,RPl,RP2 

COMMON/ CONT OR/NMINr NMAX^MMINr MMAX 

COMMON/CHR$ET/DARY,CS 

COMMON/INTSET/IrJ^NO/'MO/’DN^DM^PD 

♦initiali zati ons — — — — - 

DATA NC0L^MR0W,MARKX^MARKY/1 28^256,5,4/ 

DATA DN*DM/5*4/ 

DATA <CS(I)*I=0,1 6)/ 

' " ']• - Ml .»/.! - K« -• A« -• 

/'*'/' J"f 

♦ character array c le a i — 

DC 20 IR=OrMROW 
DO 10 IC=Or NCOL 

DARY(IC,IR)=CS(15) 

10 CONTINUE 

20 CONTINUE 

♦ character array fill 

DO 40 JaflMINrMMAX 
M0=(J“MMIN)*DM 
DO 30 I = NMIN#-NMAX 
N0=<1“NHIN)*DN 
IF(P(I/J) -LT. 3JTHEN 
PD= (“P < 1/ J ) ) 

CALL CHRFIL 
ELSE 

I F ( I -GT. NMIN)THEN 
IFCP(I-LrJ) -LT. 3) THEN 
PD = f 

CALL CHRFIL 
ENDIF 
ENDIF 

I F ( J .GT. MMIN)THEN 
IF<P(I^J-L) -LT. 3)THEN 
PD=2 

CALL CHRFIL 
ENDIF 
ENDIF 
PD=0 

CALL CHRFIL 
ENDIF 

30 CONTINUE 



40 CONTINUE 

★ printing character array ‘ 

DO 50 IR=C (MMAX-MMIN)*DM) ,0r~L 

PRINT* (132A1 ) V (OARY(IC,IRD IC=0* ( (NMAX-NMIN) *DN> ) 

50 CONTINUE 

RETURN 
END 

************************************************************************ 
SUBROUTINE CHRFIL 

★ specifications — * * - 

INTEGER CFst-XsLJ 
INTEGER N,M 

INTEGER D«DNOrMO^DN*DM*PD 

INTEGER P(Q:64*0:64),L 

REAL FS^RRMIN^RMIN^RMAX^NX^MY^MPN 

REAL F(0:64*-0:64DC(0: 1Bs1s5)#HsRf>1#RP2 

CHARACTER DARY< 0: 132#- 0: 260) *1 # C S CO: 1 6) ★ t 

COMMON FrPyC/L*H,RPlrRP2 

COMMON/ CHR SET/ DARYr CS 

COMMON/INTSET/IrJ^NOrMO/'DNrDM^PD 

★ initi all nations — * 

DATA RR MIN/ R MI N,R MAX/ 0.01*0.05^0.95/ 

★ character array D N* DM segment fill - 

LI=L 

LJ=L 

I F (PD -LT. 0) THEN 
CF=(10-PD) 

DARY(NO*HO) 3 C$(CF) 

GOTO 99 
ELSE 

I F ( PD -EQ. DTHEN 
LI=-L 
ELSE 

I F ( PD -EQ. 2) THEN 
L J=~L 
ENDIF 
ENDIF 
ENDIF 

DO 30 M^Q-r (DM-1 ) 

DO 20 N=0# ( DN-1 ) 

I F ( PD -SE. 0) THEN 

I F ( ( N+M) -LT. DTHEN 
FS=F <D J) 

ELSE 

NX=N/ REAL(DN) 

MY=M/REAL (DM) 

MPN=NX*MY 

FS=F(IxJ)*(1.0 +MPN-NX-MY)+F(I+LI/J) a(NX“MPN)+ 

: F(I, J+LJ)*(MY-MPNi+F(I+LI^J+LJ)*MPN 

ENDIF 

I F ( F S -LT- RRMINJTHEN 
C F = 1 6 
ELSE 

I F ( FS -LT. RMIN -OR. FS -GT. 

C F=1 0 
ELSE 


RMAXITHEN 



IF( F 
CF 
ELSE 
CF 
END I 
ENDIF 
ENDIF 
ENDIF 
DARYCNQ+N 
CONTINUE 
CONTINUE 
RETURN 
END 


S -GE. 0. 999)TH£N 
=INT CFS*1 0) 

=NINT(FS*10) 

F 


*LI/'HO+M*LJ) = CS(CF) 
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TABLE I. - CHARACTERISTIC IMPEDANCE AND INDUCTANCE OF 
DOUBLE LAYER MULTICONDUCTOR LINES ("INDUCTANCE" METHOD) 


■i 

Number 
of wires 

R 1 » 
mm 

r 2> 

mm 

R l, 

mm 

r 2 > 
mm 

e, 

deg 

z, 

Q 

L, 

nH/m 

8 

25 

50 

0.5 

0.5 

■9 

78.5 

396.8 

8 

25 

50 

.5 

.5 

99 

88.7 

448.2 

mm 

25 

50 

0.5 


0 

36.2 

99 

mm. 

25 

50 

.5 

19 

25.5 

37.4 

99 

32 

25 

50 

0.5 

E9 

0 

|9 

76.9 

64 

25 

50 

.5 

Mi 

0 

ESI 

29.7 

■91 

m 

25 

19 

0.5 

0 

65.6 

331 .9 

mm 

o 

25 

19 

.5 

0 

28.6 

144.4 

16 

25 

El 

0.5 

1 .0 

0 

9IM 

99 

32 

25 

99 

.5 

2.0 

0 

■■MB 

99 

32 

25 

El 

1 .0 

2.0 

0 

99 

51 .9 

64 

25 

EE 

.5 

2.0 

0 

■ESI 

20.8 

16 

10 

25 

0.5 

19 

0 

19 

87.2 

32 

10 

25 

1 .0 

m 

0 

89 

21 .8 


TABLE II. - Z AND L OF DOUBLE RING MULTICONDUCTOR 
o 

LINES ("CAPACITANCE" METHOD) 


R 1 > 
mm 

r 2> 

mm 

e, 

deg 

n 

Total number 
of wires 
in line, 

N 

Zo, 

Q 

L, 

nH/m 

25 

1 

■ 

0 

0 

.5 

8 

73.5 

371 .5 

25 



45 



8 

83.7 

423.0 

25 

I 

1 

0 



16 

33.7 

170.0 

25 

1 

a 

22.5 


’ 

16 

34.9 

176.5 

25 



■ 

0 

.5 

32 

13.9 

70.5 

25 

■ 

m 

m 

1 

1 


32 

10.4 

52.6 

25 


■ 

9 

■ 

2 


32 

6.6 

33.3 

25 

1 

1 

9 

9 

4 


32 

0.22 

1.09 

25 

50 

9 


1 


64 

3.3 

16.8 

25 

50 

9 

■ 

2 


64 

.17 

.88 

10 

25 


■ 

1 


16 

18.7 

94.4 

10 

25 

9 



.5 

32 

9.6 

48.6 

91 

25 

0 


1 


32 

■9 

27.6 

91 

25 

0 



.5 

64 

89 

15.1 



























































































TABLE III. - 64-WIRE, 5-RING TRANSMISSION LINES 


Ring 

Number 

B 

Coordinates (1st) 

Polarity 

Ring 


in ring 




(1st) 

polarity 



HH 

R, 

r, 





II 

mm 

deg 




Line 1 (all radii the same #18 AWG) 


B 

8 

0. 

51 

2.6393 

0 

+ 

A1 

t. 


8 



4.6593 

■ ■ 


- 



■9 

16 



6.6793 



+ 



n 

16 



8.6993 


■ 

- 



SK 

16 



10.7193 

B 

■ 



’ 


Line 2 (five different wire sizes) 


1 

8 

0.3 

2.6393 

0 

+ 

Alt. 

2 

8 

.51 

4.6593 

0 

- 

Alt. 

3a 

4 

.27 

6.6793 

0 

-t- 

Same 

3b 

8 

.45 

I 

22.5 

- 

Same 

3c 

4 

.51 

T 

45 

+ 

Same 

4a 

4 

.45 

8.6993 

0 

- 

Same 

4b 

8 

.35 

| 

22.5 

+ 

Same 

4c 

4 

.3 

T 

45 

- 

Same 

5 

16 

.51 

10.7193 

0 


Alt. 


DATA OUTPUT 


Line 1 

Line 2 

L, nll/m 

15.33 

18.70 

C, nF/m 

1.67 

1.37 

2 , Q 

3.03 

3.70 

A, mm 2 

26.15 

18.95 


26 



































TABLE IV. - RECOMMENDED WIRE ARRANGEMENTS FOR POWER LINE 


Number 

Wi re 
size 
(#AWG) 

Total 
number 
of wires 

Ranges of values 
inductance, 
nH/m 

Characteristic 

impedance, 

ft 

— 

Resist- 

ance, 

ft 

Weight 

of 

copper, 

kg 

High 

Low 

High 

Low 

1 

18 

80 

6.3 

10.4 

1.2 

2.1 

0.227 

87.8 

2 

18 

52 

9.7 

16.0 

1.8 

3.2 

.349 

57.1 

3 

21 

104 

4.8 

8.0 


1.6 

.342 

57.0 

4 

22 

132 

3.8 

6.3 


1.3 

.348 

56.8 

5 

24 

208 

2.4 

4.0 

0.5 

0.8 

.341 

56.8 

6 

24 

256 

2.0 

3.2 

0.4 

0.7 

.277 

69.9 


TABLE V. - DC AND AC RESISTANCE AND WEIGHT 
FOR MINIMUM NUMBER OF WIRES PER LINE 


[Line Length, 150 m (total length for resistance 
calculation is 300 m) , i.e., wires/direction.] 


Number 

AWG 

Minimum 
number of 
wi res/line 

R dc» 

ft 

^ac * ** 
ft 

Rp> 

ft* 


18 

26 

0.242 


I 

57.1 

21 

52 

.242 

mm 

HI 

57.0 

22 

65 

.244 

.245 

.343 

56.4 

33 

833 

.244 

.244 

.342 

56.4 

36 

1670 

.244 

.244 

.342 

56.4 


*Rp includes proximity effects. 

**Weight (copper alone). 
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Figure 4. - Parameter definition for equation (5). 
The direction of currents in the primed con- 
ductors is opposite to that in the unprimed. 
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Figure 5. - 64-Wire, 5-ring line 2 (inductance values in nH/m). 
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Figure 6. - 80-Wire (4x20) configuration. (All dimensions in inches.) For £ r =2.3, L = 10. 42 nH/m, C=2450pF/m, and Z = 2. 06 n. 
Scale: 0.2 inch = 1mm. 
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Conditions 

Line 

X l> 

mm 

yi* 

mm 

Polarity 

Number of 
conditions 

Spacing, 

mm 

Line 

polarity 

1 to 20 

1 

0 

0 

+ 

20 



Alt. 

21 to 40 

2 



2 








41 to 60 

3 



4 

+ 







61 to 80 

4 



6 









y = 6 ©) (m) (m) ( 204 ) ( 227 ) ( 209 ) ( 221 ) 

-61 + 

y = 4 ©) (m) © (S) ©) ( 202 ) ©) 

-41 - + 

y = 2 (l%) (m) ©) ©) (© (© ©) 

-21 + + - + 
y = 0 ( 310 ) (© (© ( 201 ) ( 227 ) (© (S) 

+1 - + - + - + 

x =0 





-20 


x = 38 


Figure 7. - 80-Wire (4x20) configuration with listing of individual self inductances in nH/m ("Inductance" method). 
L = 10.4196 nH/m ; C = Z4526 nF/m ; Z = 2. 0611 Q ; A = 32.6851 mm’ r =0.51 mm (#18 AWG). 
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Figure 8. - 80-Wire line ("Capacitance method") individual inductances. L - 6.3115 nH/m : C * 4.0490 nF/m ; 
Z 0 - 1. 2485 A - 32. 68 mm 2 ; r - 0. 51 mm (#18 AWG). 


31 




. 

Actual size 

V / 


Condition 

X l, 

mm 

n> 

mm 

Spacing, 

mm 

Polarity 

Line 

polarity 

Radii, 

mm 

1 to 32 

0 

0 

1.0 

+ 

Alt. 

0.25 

33 to 64 



1.0 



- 





65 to 96 



2.0 



+ 





97 to 128 



3.0 



- 





129 to 160 



4.0 



+ 





161 to 192 



5.0 



- 





193 to 224 



6.0 



+ 





225 to 256 


LJ 

7.0 










-33 -64 -35 +36 - + + - + - + - + - + -63 464 



+1 -2 43 -4 45 -6 4-7 -8 49 -lo +11 -12 +13 -14 +15 -16 +31 -32 


Figure 9. - 256-wire (8x32) arrangement with individual self inductances in nH/m. L = 3.2193 nH/m; C = 7.9382 nF/m; Z * 0.6368 Q ; A ■ 25. 1327 mm^ 
r - 0.25 mm ; £ r = 2.3. 
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(T) Extruded cross linked 
polyethylene multiwire 
cable (XPE) 

( ? ) Cable sealing ring 


r inner sealing 
! gasket 


r Compressible sealing 



(D Shallow trough connecting all holes, used 
for injection of sealant (Teflon or XPE) 

3B) Vacuum trough connecting all holes 

© Back-connector manufactured in three 
main pieces, molded together here 

© Crimpable pin 

© Spring tension socket (Cd HD Au or 
nickel plate) 

© 30 Percent GF, 25 percent GF glass-filled 
peek, ETFE, or similar fire rated polymer 


v 20-Gauge pins 
connected to 
female connector 


Stripped wires 
fed through to 
here and cut 


Figure 10. - Male-female junction boxes with conventional outlets. 
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30° Isometric view 

(a) High frequency power distribution ’ouss. Top cover not shown; 
material: chassis 0.063 Al-Alloy 5053; scale, 112 . 



Section A -A 

(b) Junction-box detail. (All dimensions in inches. ) 
Figure 11. - junction box. 
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